Kinetic instabilities that limit /3 in the edge of a tokamak plasma: a picture 

of an H-mode pedestal 



D Dickinsoni'^ C M Roach\ S Saarelnia\ R Scannell\ A Kirk^ and H R Wilson^ 
^ EURATOM/CCFE Fusion Association, Culham Science Centre, Abingdon, Oxfordshire, OXI4 3DB, UK 
^ York Plasma Institute, Dep't of Physics, University of York, Heslmgton, York, YOlO 5DD, UK 

March 16, 



(Dated: 



2012) 



Plasma equilibria reconstructed from the Mega-Amp Spherical Tokamak (MAST) have suf- 
cient resolution to capture plasma evolution during the short period between edge-localized modes 
(ELMs). Immediately after the ELM steep gradients in pressure, P, and density, rie, form pedestals 
close to the separatrix, and they then expand into the core. Local gyrokinetic analysis over the 
ELM cycle reveals the dominant microinstabiUties at perpendicular wavelengths of the order of the 
ion Larmor radius. These are kinetic ballooning modes (KBMs) in the pedestal and microtearing 
modes (MTMs) in the core close to the pedestal top. The evolving growth rate spectra, supported 
by gyrokinetic analysis using artificial local equilibrium scans, suggest a new physical picture for 
the formation and arrest of this pedestal. 
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Introduction:- In the high confinement mode (H- 
mode) of operation in tokamaks [l[, containment 
of magnetised plasma is enhanced through the 
formation of a narrow insulating layer that sup- 
ports steep gradients of density rie and/or temper- 
ature, at the edge of the plasma. A steep pressure 
pedestal then grows until it either saturates, or 
it is rapidly destroyed by an edge-localised-mode 
(ELM) [Ij. Pedestal growth and destruction by 
the ELM usually occurs as a cyclical process in H- 
mode plasmas. Pedestal properties are important 
in determining the capacity for plasma devices to 
confine energy: e.g. predictions of fusion power in 
the next step tokamak ITER, based on models of 
turbulent transport in the core, are sensitive to the 
size of the edge temperature pedestal 0, 0] . 

Analysis of data from many tokamaks, often 
guided by ideal MHD at infinite toroidal mode 
number, n, is consistent with the hypothesis that 
kinetic ballooning modes (KBMs) limit the pres- 
sure gradient, dP/dr, in the pedestal This 
is a key ingredient in the pedestal model called 
EPED |9l, which proposes: drift wave turbulence is 
suppressed in the pedestal by sheared flow; KBMs 
constrain the pedestal dP/dr ; and the pedestal 
width broadens as the edge current density rises, 
to trigger a finite n ideal MHD peeling-ballooning 
mode corresponding to the ELM [l3|- EPED pre- 
dictions of pressure pedestal height and width prior 
to the ELM crash, are consistent with data from 
several tokamaks 0, Q , and partially supported by 
recent data from the spherical tokamak MAST 
[llj . The EPED model, however, cannot fully ex- 
plain the evolution of kinetic profiles. 

This Letter presents detailed kinetic analyses of 
an evolving pedestal equilibrium from MAST Q, 
and probes whether microinstability physics helps 
explain: (i) the inwards expansion of steep pedestal 
profile gradients into the core, and (ii) mechanisms 
that limit the fully developed pedestal. 
Inter-ELM Pedestal Profile Evolution in MAST:- 



Data were taken from reproducible MAST H-mode 
discharges with regular type I ELMs, and heated 
by 3.4MW of neutral beam injection Q. Profiles 
of electron density, rie, and electron temperature, 
Te, were measured every 4.2ms using a Thomson 
scattering (TS) system with 130 spatial points [12 1. 
Spatial resolution was sufficient to resolve the steep 
pedestal profiles on the high field side. Figure [TJa- 
c) shows TS profiles from one ELM cycle that 
demonstrate stronger changes in than T^, with 
the density pedestal growing throughout the cycle. 
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FIG. 1. TS profiles of (a) Te and (b) rie from three 
times in one inter-ELM period that are indicated on 
the D-alpha trace (c). TS data from periodic ELMs 
were used to obtain (d) Pe{^ N,t) over the ELM cycle. 

Time resolution of the ELM cycle was enhanced 
by ordering TS profiles, from the periodic ELM 
phases of the discharges, in time with respect to the 
previous ELM. This revealed that dPe/dr recovers 
very rapidly close to the separatrix, during the first 
10% of the cycle after the ELM crash. These pro- 
files were used, assuming ion temperature, Ti = Te, 
to reconstruct equilibria (see Q for full details), ac- 
commodating the edge bootstrap current [l3[, at 
five times normalised to the inter-ELM period, t, 
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during one ELM cycle. Figure [Ijd) shows how the 
edge electron pressure, Pe(^'Af), evolves (where \I'Ar 
is normalised poloidal flux): changes in Ue domi- 
nate Pe evolution; the Pe pedestal top moves in- 
wards from ^ 0.98 at t = 0.1 to ~ 0.96 at 
t = 0.9; and the edge region of high dPe/dr (i.e. 
the edge transport barrier) expands into the core 
with modest change in the peak dPe/dr. Similar 
observations during regular type I ELMs have been 
reported by the tokamak experiments DIII-D 14] 
and ASDEX Upgrade 

Micro stability during Pedestal Evolution:- The H- 
mode pedestal is associated with a strong pressure 
gradient, which can drive electromagnetic instabil- 
ities, like the KBM, where the magnetic perturba- 
tions are crucial to the instability mechanism. Col- 
lisionless KBMs and ideal MHD infinite-n balloon- 
ing modes are described by related equations and 
have similar character Il6j . Finite Larmor radius 
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(FLR) effects, only included in the kinetic treat- 
ment, are often stabilising [l7| : the region unsta- 
ble to infinite-n ideal MHD ballooning modes can 
be broader than the region unstable to KBMs [l3| . 
The stability of infinite-n ballooning modes is of- 
ten, nevertheless, found to be a reasonable proxy 
for KBM stability Q . The gyrokinetic equation is 
derived from a first order expansion of the Vlasov 
equation in the small parameter S = p/L (where 
X is a typical equilibrium gradient scalelength) . 
Short equilibrium gradient scalelengths increase S 
in the pedestal: for ions Si ~ 0.3 in the steep- 
est part of the MAST pedestal, but Si is an order 
of magnitude smaller immediately inside the fully 
developed pedestal top and for electrons Se ^ 1 
throughout the plasma. Ideal MHD is applied rou- 
tinely to analyse pedestal equilibria. While gy- 
rokinetics is less accurate in the pedestal than in 
the core, it improves on MHD in the high n limit 
through the inclusion of FLR corrections and reso- 
nances. Gyrokinetics was previously used to study 
conditions where electron drift waves may cause 
transport in a siniple slab geometry model of the 
H-mode pedestal flSj . 

Local gyrokinetic analysis, using the initial value 
fiux-tube code GS2 [lH , has obtained the microsta- 
bility properties of the MAST edge plasma during 
the ELM cycle. GS2 solves the electromagnetic 
gyrokinetic equation [2^ for each plasma species, 
to find the fastest growing microinstability and its 
growth rate, 7, for specified values of kyPi (where 
ky is the in-flux-surface component of the perpedic- 
ular wavenumber). Dominant growth rate spectra 
covering 0.05 < kyPi < 5.5 were computed on 12 
surfaces spanning 0.94 < '^n < 0.995, for 5 times 
during the ELM cycle, t = (0.1,0.3,0.5,0.7,0.9). 
Sheared toroidal fiow was neglected in this first 
analysis, and is not expected to be important just 
inside the pedestal top. 

In the steep dP/dr pedestal region, the domi- 
nant microinstabilities grow around kyPi ^ 0.2 and 
are KBMs with twisting parity. There is a sharp 



FIG. 2. Real part of the dominant eigenfunctions of Ay 
at t = 0.5 plotted as functions of ballooning angle, 0, 
exhibiting (a) tearing parity at kyPi — 3.28, ^'jv ~ 0.95 
and (b) twisting parity at kypi = 0.218, ^^jv = 0.98. 

transition to microtearing modes (MTMs) growing 
around kyPi ~ 3 in the shallower gradient region 
immediately inside the top of the pedestal. The 
evolution of the radial region where KBMs dom- 
inate is similar to the region found to be unsta- 
ble to infinite-n ideal MHD ballooning modes Q. 
Figure [2] illustrates typical eigenfunctions of the 
parallel magnetic vector potential, for MTMs 
and KBMs. KBMs are driven by both tempera- 
ture and density gradients, propagate in the ion 
diamagnetic drift direction, and typically have a 
ratio of electron collision frequency to mode fre- 
quency z^e/w > 1. MTMs are driven by the elec- 
tron temperature gradient at large r]e = Ln^/Lx^, 
propagate in the electron diamagnetic drift direc- 
tion, and have electron collisionality z^e/w ^ 0(1). 
The MTMs have similar characteristics to modes 
reported at mid-radius in spherical tokamak plas- 
mas j2T| - [23j . with two notable exceptions: the edge 
modes here (i) are considerably less extended in 
ballooning angle, 0, and (ii) arise at substantially 
higher kyPi. Microtearing modes have also pre- 
viously been reported to be unstable close to the 
edge in ASDEX Upgrade and in conceptual 
spherical tokamak burning plasma devices [25l . [26j . 
Recent measurements from the pedestal region be- 
tween ELMs in DIH-D [13 find high and low uj 
bands of turbulence propagating in the electron 
and ion drift directions respectively, which are con- 
sistent with the properties of the MTMs and KBMs 
reported here. 
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FIG. 3. ^{kypi,t) over the ELM cycle at *]v = 
0.95,0.96,0.97,0.98 in (a)-(d) respectively. 7 is nor- 
malised to Vti/a, vti is the local ion thermal velocity 
and a is the minor radius. 

Figure [3] shows the ELM cycle evolution of 
the growth rate spectrum on four flux surfaces 
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*Ar = 0.95,0.96,0.97,0.98. Deepest inside the 
core, ^'at = 0.95 sits just inside the shallow gra- 
dient core plasma when the pedestal is fully de- 
veloped (corresponding to surface b of Figure S] 
below). MTMs dominate the spectrum in Fig- 
ure Elja) , which broadens and increases in ampli- 
tude as locally /3 rises through the MAST ELM 
cycle. MTMs may limit the electron tempera- 
ture gradient in this region of the plasma. We 
note that magnetic drifts are important to the 
microtearing mode [2l| . and that dT^/dr may be 
more severely limited in equilibria with more un- 
favourable drifts. Further out, Figure [3l[c) shows a 
striking transition in j{kyPi), when the '^n = 0.97 
surface joins the expanding pedestal at t ^ 0.6: 
the maximum growth rate, -^"1^''^ falls and then 
rises; and the dominant wavcnumber bifurcates to 
a lower value. Inspection of the eigenfunctions re- 
veals that MTMs, centred at kyPi ^ 3, dominate 
while ^'at = 0.97 experiences core conditions, but 
that the MTMs are suppressed and supplanted by 
KBMs at kyPi ~ 0.2 when dP/dr steepens to the 
higher pedestal value. The physics driving this 
transition is investigated below. Further into the 
pedestal at ^'at = 0.98, Figure |3|d) shows that 
KBMs are driven unstable earlier, but their growth 
rates fall later in the ELM cycle. KBMs remain the 
dominant modes, but for t > 0.5 increasing boot- 
strap current reduces the magnetic shear, which 
pushes the modes towards marginal stability. 

How might the microinstabilities physically in- 
fluence the evolution of the MAST pedestal? To 
deepen our understanding of the basic mecha- 
nisms, several artificial stability scans have been 
performed based around the key surfaces illus- 
trated schematically in figure |H 
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FIG. 4. Schematic P profiles from during the growth 
of the pedestal (dashed green) and when it is fully de- 
veloped (solid red), indicating three reference surfaces 
for stabihty scans. 

Why does the pedestal expand inwards?- Artifi- 
cial scans, around the reference equilibrium flux- 
surface a in Figure HI have probed in more detail 
the mechanism by which the pedestal expands in- 
wards. Surface a (^'at = 0.97,t = 0.5) sits just in- 
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FIG. 5. 7(fc yPi) scan around surface a {'^n = 0.97 
t — 0.5) scaling dP/dr and R/Ln self-consistently. 



side the shallow gradient core plasma, and will im- 
minently be subsumed by the advancing pedestal. 

Figure [5] shows the growth rate spectrum from 
an equilibrium scan around surface a, multiplying 
the pressure gradient /?' by a factor (3'^^^ and scaling 
the density gradient R/Ln and 77e self-consistently 
at fixed Tie, Te, Ti, and temperature gradients 
R/Lt^ and R/Lt^- This scan closely resembles the 
evolution observed during the MAST ELM cycle, 
where the steepening pressure gradient is domi- 
nated by an increasing density gradient. The dom- 
inant 7 and ky fall as /3' increases above the exper- 
imental value. The tearing parity mode is predom- 
inantly stabilised by the disruption of the relative 
phases between current and density perturbations 
arising from the increasing R/Ln drive terms, and 
further stabilisation arises from more favourable 
drifts at higher j5' [H, [2^. This stabilisation of 
dominant MTMs reduces all transport associated 
with these modes. The quasilinear electron heat 
flux from the modes is larger than the particle flux, 
as is typical for MTMs. At higher R/Ln, KBMs 
become dominant as their growth rates increase; 
KBMs transport heat and particles in all channels. 
The maximum /3' at ^fjv = 0.97, from late in the 
ELM cycle, corresponds to jS'^^^ = 1.78 which lies 
in the KBM region of figure [5j 

This scan suggests that: electron heat flux, 
and modest particle flux, initially fall as drig/dr 
increases, due to MTMs becoming suppressed; 
and at higher drie/dr all transport fluxes increase 
strongly when KBMs become dominant. This re- 
duction in transport will allow the electron pres- 
sure and density pedestals to build up during the 
ELM cycle. The particle source from edge neutrals 
is likely to play a significant role in this evolution of 
Ue ■ This picture is qualitatively consistent with the 
inwards advance of the pedestal in MAST, whereby 
the shallow pressure gradient at a in Figure S] un- 
dergoes a rapid transition to the steeper value in 
the pedestal. The relative impact on the evolution 
of He and Tg pedestal profiles will be sensitive to 
edge sources of heat and particles, and most no- 
tably the particle source from cold neutrals. A 
transient increase in the particle source, e.g. from 
a pellet, may accelerate the development of the 
density and pressure pedestal. 

What limits the inward expansion?:- Surface b of 
Figure IK^-jv = 0.95, t = 0.9) sits wcU within the 
shallow gradient core plasma when the pedestal is 
almost fully developed. (3 would clearly have to 
increase on this surface if the pedestal were to ad- 
vance further into the core by raising dPe/dr on 
surfaces where 'I'Af > 0.95. This motivated an ar- 
tificial scan centered on surface b, where /3 and rig 
were scaled by /3fac consistently with the observed 
pedestal evolution in MAST. /?' (which influences 
the magnetic drifts) was adjusted self-consistently 
for fixed dn^/dr, T, R/Lt and collision frequen- 
cies. The width of the unstable MTM spectrum 
was found to increase marginally with /3fac. 
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FIG. 6. 7 spectrum scan around ^^jv = 0.96, t = 0.9 
(surface c) scaling /3 and rie as described in the text. 

More interestingly, this scan was repeated on 
surface c of Figure |4]( 5' AT = 0.96, t = 0.9), which is 
further out than &, in the transition region between 
the shallow core and the steep pedestal pressure 
gradient just prior to the ELM. Figure H] indicates 
that this surface sits close to a critical value of 
/3 where strongly growing KBMs and MTMs 
coexist over a broad spectrum and are equally 
dominant at the arrow, where n ~ 25. The cross- 
ing of this threshold may trigger a significant trans- 
port event in the edge plasma, n = 25 coincides 
with the toroidal mode number of the peeling- 
ballooning mode found marginally unstable at the 
end of this ELM cycle Q, and is in the observed 
range for ELM filaments in MAST H^. 
Conclusions:- Plasma equilibria have been recon- 
structed from the spherical tokamak MAST, with 
sufhcient resolution in time and space to capture 
plasma evolution during the short period between 
ELMs. Immediately after the ELM, steep pedestal 
pressure and density gradients form close to the 
separatrix, and these advance into the core. Fully 
electromagnetic local linear gyrokinetic analysis re- 
veals the dominant microinstabilities 0jS &- function 



of radius and time through the ELM cycle. This 
suggests a new physical picture of the formation of 
this pedestal: steep pedestal pressure and density 
gradients are limited by KBMs; MTMs dominate 
and limit R/Lt^ in the shallower pressure gradi- 
ent region of the core just inside the pedestal; the 
pressure pedestal can propagate into the core be- 
cause increasing R/Ln and dP/dr stabilizes the 
MTMs until they become supplanted by KBMs at 
a higher pressure gradient; and meanwhile in the 
core MTMs become increasingly virulent as /? in- 
creases in the shallower pressure gradient region. 
Furthermore, when the pedestal is nearest to be- 
ing fully developed, the gradient transition region 
is close to an instability threshold where MTMs 
and KBMs become simultaneously unstable with 
large growth rates over a broad spectral range. 
Breaching this limit may trigger signicant change 
in the edge transport. This picture should be com- 
pared with data from a broader range of H-mode 
plasmas and from turbulence diagnostics, and ex- 
tended to predict transport timescales using non- 
linear gyrokinetic fluxes and realistic sources. 
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